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We have developed a unique micromechanical method to extract extremely thin graphite samples. 
Graphite crystallites with thicknesses ranging from 10 - 100 nm and lateral size ~ 2 /im are extracted 
from bulk. Mesoscopic graphite devices are fabricated from these samples for electric field dependent 
conductance measurements. Strong conductance modulation as a function of gate voltage is observed 
in the thinner crystallite devices. The temperature dependent resistivity measurements show more 
boundary scattering contribution in the thinner graphite samples. 



The recent focus of research in graphitic materials, 
such as fullerenes 1} and carbon nanotubes [J 0] has 
renewed interest in electronic transport in graphene, a 
single atomic layer of graphite. While graphite is a 
semimetal with strong electron and hole compensation, 
graphene is expected to be a zero-gap semiconductor with 
a vanishing density of states at the charge neutral point. 
Electrical transport in thin graphite crystals composed of 
only a few graphene layers is of particular interest in elu- 
cidating the evolution of electronic structure from bulk 
single crystals to two-dimensional planar systems. Ex- 
perimental work to synthesize very thin graphitic layers 
directly on top of a substrate or to extract using chem- 
ical Q , or mechanical 0, @ pathways has been demon- 
strated to produce graphitic samples with thicknesses 
ranging from 1 - 100 nm. Systematic transport mea- 
surements have been carried out on mesoscopic graphitic 
disks and cleaved bulk crystals 8] with sample thick- 
nesses approaching ~ 20 nm, exhibiting persisting bulk 
graphite properties at these length scales. However, elec- 
tric field dependent charge transport has not been mea- 
sured, due to the fact that the size of the samples and 
the geometry of the devices are inadequate for such a 
measurement. 

In this letter, we present a unique method to extract 
mesoscopic graphite crystallites, with thicknesses rang- 
ing from 10 - 100 nm, from bulk highly oriented py- 
rolytic graphite (HOPG) using micromechanical manip- 
ulation and microfabrication techniques. Mesoscopic de- 
vices consisting of these graphite crystallites are fabri- 
cated for electrical transport measurement. A degener- 
ately doped Si substrate serves as a back gate, allow- 
ing us to probe the electrical properties as a function 
of charge concentration. Electrical conductance of these 
mesoscopic graphite devices exhibits significant gate volt- 
age dependance. From the temperature dependent re- 
sistivity measurements in these samples, we deduce the 
scattering rate of mesoscopic graphite of varying thick- 
nesses. 

The graphite samples used in this study are extracted 
from HOPG (Grade-1, Structure Probe, Inc.). Follow- 
ing a similar technique demonstrated in [(|, arrays of 
graphite micro-pillars (~ 2 x 2 x 5 /im 3 ) are fabricated 
on the HOPG surface using micropatterning followed by 
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FIG. 1: (a) scanning electron microscope image of an HOPG 
crystallite mounted on a microcantilever. Inset: Bulk HOPG 
surface patterned by masked anisotropic oxygen plasma etch- 
ing, (b) Schematic drawing of the micro-cleaving process, (c) 
Thin graphite samples cleaved onto the Si02/Si substrate, 
(d) A typical mesoscopic device fabricated from a cleaved 
graphite sample. 



masked anisotropic oxygen plasma etching (Fig. ^a) in- 
set). Once the array is formed, an individual block of 
the pillar is removed from the surface using a precision 
micro-manipulator under a high resolution optical micro- 
scope. The detached graphite block is then transferred to 
a micro- machined silicon cantilever (MikroMasch) , where 
it is glued down by a small amount of ultra violet sensi- 
tive epoxy (Fig. IHa)). 

In the next step, we use the mounted graphite block 
on the cantilever as the tip of an atomic force micro- 
scope (AFM) in order to transfer thin graphite samples 
onto a Si02/Si substrate for subsequent device fabrica- 
tion. By operating the AFM (Nanoscope Ilia Multimode, 
Digital Instruments) in contact mode with load on the 
graphite mounted cantilever, very thin layers of HOPG 
are sheared off onto the substrate (Fig^b)). This mi- 
croscopic cleaving process can be controlled by tuning 
the loading normal force between the cantilever and the 
substrate. The van der Waals binding energy between 
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FIG. 2: Normalized conductance measured in samples of vary- 
ing thickness as a function of gate voltage at T =1.7 K. Inset: 
Calculated density of states for bulk graphite [TfT ]. 

the graphene layers in graphite is ~ 2 eV nm -2 As- 
suming the graphite-substrate friction coefficient is ~ 1, 
the required normal force to cleave off a 1 /zm 2 graphite 
sheet from the top of the crystallite is ~ 300 nN. A nor- 
mal force in the range of 10 - 2000 nN can be easily 
managed by commercially available silicon cantilevers. In 
our experiment, by fine tuning the normal force, square 
graphite crystallites with lateral size ~ 2 /jm and thick- 
nesses ranging from 10 - 100 nm are easily obtained 
using this process (FigGfc)). Once the graphite sam- 
ples are cleaved onto a substrate with predefined align- 
ment marks, AFM images are acquired for measuring the 
height of the cleaved crystallites and registering the po- 
sition for device design. The Cr/Au (1 nm/30 - 180 nm) 
electrodes are then fabricated by electron beam lithogra- 
phy. Typically multiple contacts are fabricated around 
the square graphite crystallites allowing us to perform 
conductivity measurements using van der Pauw method. 
In addition, a degenerately doped silicon substrate is used 
as a gate electrode with 500 nm of thermally grown sili- 
con oxide acting as the gate dielectric. 

Fig. |21 shows the conductance, G, versus the gate 
voltage, V g , at 1.7 K for samples with varying thick- 
nesses. The curves are normalized with their minimum 
values G min = G(V g min ), where V™ in « -40 V. As 
AV g = Vg—V™™ moves away from zero, G(V g ) increases, 
increasing more rapidly for thinner samples. Specifically, 
at AV g = 100 V, G/G min is 1.47 and 1.12 for the samples 
with thickness 12 and 42 nm respectively, while no ap- 
preciable change is observed in the 95 nm sample. This 
observation implies a larger electric field effect in thinner 
samples, as expected. 

Considering the electrostatic coupling between the 
samples and the gate, the observed behavior of G(V g ) 
at low temperatures can be understood quantitatively as 
discussed below. When the gate voltage is applied, the 



electrostatic potential at the surface of the sample rises to 
00) and a charge density rii n d(z) is induced in the sample, 
where z is measured from the interface of the sample and 
the substrate. Following the Thomas-Fermi approxima- 
tion |ll| , the electrostatic potential in the sample is given 
by 4>[z) — 4>Qe~ z / Xs , where 4>$ is a constant (determined 
below) and A s is the screening length. \ s can be inferred 
from 0/A 2 = —eriind/ 'eo, where e and £o are the electron 
charge and vacuum permittivity, respectively. The total 
amount of induced charge in the sample with thickness 
d is related to the capacitance of the sample to gate per 
unit area, C g , by C g (V g - <fi ) = - en ind (z)dz. Eval- 
uating this integral, we obtain </> = V g /(1 + Cg/Cg), 
where the sample dependent quantum capacitance per 
unit area Cq — sq(1 — e~ d / Xs )/ \ s . 

In order to estimate the enhanced conductance in 
the presence of an induced charge density gradient, we 
use Einstein's relation for local conductivity <r{z) = 
e 2 N(E F )D(E F ), where N(E F ) and D(E F ) are the den- 
sity of states and and the electron diffusion coefficient 
at Fermi energy E F , respectively. In graphite both of 
these quantities are sensitive to the change of E F , thus 
it should be accounted for in the determination of lo- 
cal conductivity in the presence of an electric field. As- 
suming our samples retain the band structure of bulk 
graphite, the functional form of N(e) can be obtained 
from graphite band structure calculations (Fig. |3 in- 
set) [lfjj. Here e is measured from the charge neu- 
tral point, where N(e) has its minimum value. For 
the energy range |e| < 300 meV, we can approximate 
N(e) w N (l+j3\e\), where N = 5.2xl0 20 cm-^V 1 
and {3 = 4.8 ev _1 . A further assumption is required 
to estimate D(E F ) in the presence of an electrostatic 
potential. Here, we adopt the simple two band (STB) 
model [T3 , which has been successful in explaining elec- 
tron transport in bulk graphite |13| . Assuming the scat- 
tering time, r, and effective mass of carriers are rela- 
tively insensitive to change of E F , the STB model yields 
D(e) cx v f t oc \e + Eq\, where v F is the Fermi velocity 
and Eo half the band overlap between the electron and 
hole bands. With these assumptions, the normalized con- 
ductance of the sample is given by 

G(V g )/G mm = I dzN{e)\e + E \/N Q E d (1) 
Jo 

where e(z) — —e4>oe~ z / Xa . In order to compare this equa- 
tion to our experimental observations, the values of X s 
and E are needed. We use the following values from 
the literature: X s = 0.4 nm 14] and E = 0.015 eV pj . 
Employing these values, we estimate C g by comparing 
Eq. ^to the experimentally observed values. The corre- 
sponding C g are 30 aF/zm -2 and 27 aF/im~ 2 for 12 nm 
and 42 nm thick graphite devices, respectively. These 
values are in reasonable agreement with the value of C g 
obtained by considering the device geometry. 

We now turn our attention to the temperature depen- 
dence of the electrical conductance in samples of vary- 
ing thickness. Fig. [3] shows the resistivity normalized to 



3 




100 200 300 

T[K] 

FIG. 3: Resistivity (normalized to its 300 K value) as a func- 
tion of temperature for different thicknesses. Broken lines are 
fits with the STB model. Inset: The mobility due to static 
disorder, as a result of fits to the STB model. The dotted line 
is an empirical fit. 

its room temperature value, p(T) / p(300K ), as a func- 
tion of temperature, T, for d ranging from 12 - 95 nm at 
V g = 0. While the thicker samples show a metallic behav- 
ior (dpjdT > 0), thinner samples develop a more com- 
plicated behavior in T . This behavior can be understood 
by considering two competing factors as the temperature 
varies. Assuming the mobility, p, is similar for the com- 
pensating electrons and holes in the samples, the sample 
resistivity is given by p — l/(\e\pn), with the total charge 
carrier density n = n e +n/i, where n e and are the car- 
rier densities of electrons and holes, respectively. In the 
STB model, for completely compensated samples (i.e., 



AV g = 0), n/2 = n e = n h = C fc B Tln(l + exp^), 
where kg is Boltzmann constant and Co a constant inde- 
pendent of T ^3] • We further assume that the mobility 
can be expressed as p~ x = nj + AqT, where pj is 
the contribution from static scattering centers and Aq 
is a constant depending on electron-phonon scattering 
in graphite. Then, as T increases, p decreases due to 
increasing electron-phonon scattering, while n increases 
rapidly due to the thermal excitation of charge carriers 
around Ep. Using E and A pi as fitting parameters, 
this model fits well with our experimental observations 
over a large temperature range and samples with varying 
thicknesses (broken lines in Fig. [SJ. As a result of the 
fit, we obtain Eq ~ 10 - 20 meV, which is in reasonable 
agreement with the values found in the literature [15J . In 
addition, we obtain A^pi from the same fit. Assuming 
Aq is similar for all measured samples, the pi for different 
samples can be compared. The inset of Fig. [21 shows de- 
creasing pi as sample thickness is reduced. An increasing 
contribution of boundary scattering at the sample sur- 
faces can explain this behavior. We note, however, that 
the reduction of pi seems to saturate as d — > 0. In ad- 
dition, our recent galvanomagnetic measurements |16| in 
these samples showed that the quantum oscillations in 
electron transport can be readily observed at high mag- 
netic fields. Both observations suggest that the quality of 
the sample has not degraded after our micromechanical 
extraction from the bulk. 
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